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Abstract. In the context of explosion models for Type la Supernovae, we present one- and two-dimensional simulations of 
fully resolved detonation fronts in degenerate C+O White Dwarf matter including clumps of previously burned material. The 
ability of detonations to survive the passage through sheets of nuclear ashes is tested as a function of the width and composition 
of the ash region. We show that detonation fronts are quenched by microscopically thin obstacles with little sensitivity to the 
exact ash composition. Front-tracking models for detonations in macroscopic explosion simulations need to include this effect 
in order to predict the amount of unburned material in delayed detonation scenarios. 
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1. Introduction 

Driven by the ever increasing amount and qual ity of observa- 
tional data for Type la Supe rnovae (SNe la) (e.g jLeonard et all 
l2005t iMazzali et alJl2005t) and their continuing role as cos- 
mological distance indicators jRiess et alJ E"005). the theoret- 
ical efforts to understand the physics of exploding C+O White 
Dwarf stars near the Chandrasekhar mass have progressed sig- 
nificantly. We now have at our disposal three-dimensional hy- 
drodynamical models including different algorithms to track 
thermonuclea r flames, d etonatio ns, and unr esolved turbulence 
(e.g.. lGamezo et all2005URopke et alJ2005h . as well as the re- 
sults of fully resolved Rayleigh-Taylor unstable flame fronts 
showing the transiti on to isotropic turbulence on small scales 
JZingale et alJ I2005I) . Nevertheless, no consensus has been 
reached on whether pure turbulent deflagrations can explain SN 
la observations as indicated by high-resolution multi-point ig- 
nition models or whether delayed detonations of some sort are 
needed to burn more material at low densities. 

It is undisputed, however, that any detonation must be pre- 
ceeded by a phase of subsonic burning that pre-expands the 
star and, as a result of the Rayleigh-Taylor (RT) instability, 
leaves behind highly inhomogeneous regions of burned and 
unburned material. Large, merging RT mushrooms frequently 
produce large volumes of fuel fully enclosed by ash (but not 
vice versa). One of the key virtues of delayed detonations is 
their suspected ability to deplete much of the remaining un- 
burned C+O in the White Dwarf's core region. Large enough 
quantities of such material would produce low-velocity carbon 
lines in the nebular spectra of SNe la whose observational ab- 
sence can already rule out coarsely r esolved deflagration mod- 
els with single-point central ignition (Kozma et al. 2005). 



The work reported in this paper sidesteps the key issue 
for evaluating the "naturalness" of SN la scenarios that in- 
voke a detonation: under which microscopic conditions can a 
thermonuc lear detonation in degenerate C+O material be ig- 
nited (e.g. , iNiemever & Wooslevlfl997l iKhokhlov et allll997t 
iNiemeven 1999h ? We instead assume that a detonation has 
formed by some unspecified mechanism and begins to prop- 
agate both outward and inward from its point of birth. Upon 
impacting on a region of ash produced by the preceeding defla- 
gration, thermonuclear reactions cease to provide the overpres- 
sure that compensates for the dissipation of the shock front. 
After propagating through the ash and entering fresh fuel on 
the other side, the shock may either be sufficiently powerful 
to re-ignite a detonation or not, in which case fully enclosed 
pockets of fuel remain unburned. 

One aim of this work was to numerically determine the lin- 
ear size of a sheet of ashes that a C+O detonation can still 
cross and survive 1 . This was done with ID, fully reso lved sim- 
ulations using the FLASH code jFrvxell et al-l feoOO ) that has 
previously been em ployed successfully to study thermonuclear 
C+O detonations by Ti mmes et alJfcOO Ob). Only for the lowest 
fuel density in our parameter space, p = 3.2 • 10 7 g cm^ 3 , did 
we find successful detonation re-ignitions up to an ash thick- 
ness of roughly 6 detonation widths (corresponding to 8 cm). 
In all other cases, the detonation failed to re-ignite on the other 
side, regardless of the detailed composition of the burned mat- 
ter (pure 56 Ni, 28 Si, or NSE). 

In addition, we carried out a small number of 2D simula- 
tions in order to test the ability of detonations to pass through 
thin funnels of fuel separating neighoring regions of ash. Not 



1 After completion of this work we were in formed that a similar 
study had previously been performed by Livne 1 2004). 
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surprisingly, we found that a minimum width of these funnels 
exists that allows the detonation to pass if this width is close to 
the detonation thickness. 

The structure of this paper is as follows. We describe the 
numerical setup and calibration runs in Sect. 2, followed by 
a summary of the ID (Sect. 3.1) and 2D (Sect. 3.2) parame- 
ter studies. Conclusions are presented in Sect. 4. , in particu- 
lar those concerning the relevance of our results for large-scale 
simulations of detonations in exploding white dwarfs. 

2. Tests of the numerical setup 

The FLASH code is a modular, adaptive mesh, parallel simula- 
tion code capable of handling general compressible flow prob- 
lems in astrophysical enviroments. It uses an explicit, direc- 
tiona lly split version of the p iecewise-parabolic method (PPM) 
JColella & Woodward! \l 984) to solve the compressible Euler 
equations and allows for general equations of state using the 
method of lColella & Glad fl 985). An equation of state appro- 
priate for handling the interior of a C+O White Dwarf is im- 
plemented using a thermodynamically consistent table lookup 
scheme dTimmes & Swe stv 2000) . Source terms for thermonu- 
clear reactions are solved using a semi-implicit time integrator 

J hi — —j 

coupled to a sparse matrix solver (Timmes 1999). Further de- 
tails about the algorithms used in the code, the co de's structure, 
and r esults of verification tests are documented in lFrvxell et alJ 
(200(3). 

As a good compromise between speed and precision, a 1 3 
isotope a-chain plus heavy-ion reaction network was used in 
our calculations. A network of this kind (called approxl3) 
using the isotopes 4 He, 12 C, 16 0, 20 Ne, 24 Mg, 28 Si, 32 S, 36 Ar, 
40 Ca, 44 Ti, 48 Cr, 52 Fe and 56 Ni is i mple mented in FLASH an d 
thorougly tested by Timmes ( 1999) and lTimmes et alJ fcoOOa). 

In order to quantify the necessary resolution for our simu- 
lation of detonation fronts with FLASH, we carried out a set 
of ID test simulations. The boundary conditions consisted of 
a system of length / = 1024 cm with a reflecting boundary at 
x = cm and an outflow boundary condition at x = 1024 cm. 
As initial conditions we chose a composition of 50% 12 C and 
50% ls O at a temperature of 10 7 K and a material velocity of 
v x = cm s _1 . To ignite the detonation we set the inital condi- 
tions in an ignition area between x = cm and x = 25.6 cm 
to a temperature of 10 10 K and a material velocity of v x = 
10 9 cms _1 to ensure sufficient numerical diffusion for a suc- 
cessful detonation ignition. We examined the speed and the nu- 
clear energy generation rate of the detonation for densities p 
between 1 .0 • 1 7 and 1 .0 • 1 9 g crrT 3 , depending on the resolu- 
tion of the simulation. The speed of the detonation was evalu- 
ated by taking the derivative of the position for the maximum of 
the nuclear energy generation rate throughout t he system over 
time. We were able to reproduce the results of Sharpe i 19991) 
except for a difference of about 1 % that can be accounted for 
by small differences in the employed reaction networks. 

As the detonation speed is not very sensitive to the numer- 
ical resolution, we also examined the magnitude of the maxi- 
mum of the nuclear energy generation rate for effective spatial 
resolutions of 0.4 cm, 0.1 cm, 2.5 • 10~ 2 cm, and 6.25 • 10~ 3 cm. 
As can be seen in Fig. the chosen inital conditions gener- 



1 0- f 





I I I I I I I I 


1 1 | 


=6 : 




i\ 

j \ 


..... ref 


=8 - 




\ 


- - ref m 


= 10 






— - ref m 


=12 




1 \ \ 

£\\ 








it ; '■ >. 

A V 

i V V 

l l T T 1 , , , 1 


I 





2x10- 7 4x10- 7 6*10- 7 8x10- 7 

t/s 

Fig. 1. Magnitude of the maximum of the nuclear energy gen- 
eration rate for a density p = 5.0 • 10 7 g crrT 3 and different 
spatial resolutions 0.4 cm ( ref mflA = 6) , 0.1 cm ( ref max = 8), 
2.5 • 10~ 2 cm ( Tsf max = 10), and 6.25 • 10~ 3 cm ( ref mflA - = 12) 

ate a detonation with a high nuclear energy generation rate that 
quickly converges to a stationary state. One can see that in- 
creasing the resolution decreases the fluctuation amplitude of 
the maximum of the nuclear energy generation rate up to the 
point where the necessary resolution for a fully resolved deto- 
nation is reached. From this analysis we found resolutions of 
0.1 cm, 2.5 • 10 -2 cm, and 6.25 • 10" 3 cm sufficient for fuel den- 
sities of 3.2 • 10 7 g cirT 3 , 3.2 • 10 7 g cirT 3 , and 1.0 ■ 10 8 g cirT 3 , 
respectively. All further simulations were performed with these 
resolutions. 

For a density of 1.0 • 10 7 g cm 4 , we found that our initial 
conditions did not generate a stable detonation, and for higher 
densities p = 3.2- 10 8 . . . 1.0 ■ 10 9 g cm 3 , we were unable to 
reach the necessary resolution for a fully resolved detonation. 

A convenient length scale to normalize our results is given 
by the detonation width b, defined here as the distance between 
the point of maximum pressure and the point where the nu- 
clear energy generation rate has dropped to 10% of its max- 
imum (Fig. |3 . We found widths b of 1.37 cm, 1.00 cm, and 
0.52 cm for fuel densities of 3.2 • 10 7 g cm -3 , 5.0 • 10 7 g crrT 3 , 
and 1.0 • 10 8 g cm 4 . Width-normalized lengths will be denoted 
henceforth by l s = l/b(p). 

3. Interaction of detonation fronts with ash 
obstacles 

3.1. Properties of the ash obstacles 

The composition and thermal state of the nuclear ash produced 
by the deflagration front preceding the detonation depends on 
the fuel density at the time it was consumed. In particular, the 
ash composition can be expected to affect the behavior of the 
detonation front since the shock energy is dissipated by the dis- 
sociation of heavy nuclei. 

In order to study the sensitivity of shock dissipation to ash 
composition, we examined five different setups (N, S, M1-M3) 
covering a wide range of possible compositions (see TableQ. 
The ash regions were constructed by changing the composition, 
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N 


S 


M2 


r/io" k 


5.505 


5.265 


5.336 


p/10 7 g cirr 3 


1.919 


2.205 


2.116 


Aq/10 11 erg g -1 


7.860 


5.977 


6.268 



Table 3. Parameters of the ashes for p = 5.0 • 10 7 g cm 
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M3 


r/io 9 k 


6.725 


6.366 


6.357 


p/10 7 g cnr 3 


4.456 


5.043 


5.027 


A^/10 17 ergg- 1 


7.860 


5.977 


6.007 



J Table 4. Parameters of the ashes for p = 1.0 • 10 8 g cm 



Fig. 2. Temperature, pressure, and nuclear energy generation 
profile of a detonation front for p = 5.0 • 10 7 g cm -3 4.0 • 10 7 s 
after ignition. The dotted lines depict the detonation width b 
according to our definition. 
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Fig. 3. Initial conditions for a system with p = 5.0 ■ 10 7 g cm -3 
and an ash obstacle type M2 with width w = 128 cm . From top 
to bottom: temperature, density, pressure. 



Table 1. Composition of different types of ash. The ashes of 
type Ml, M2, and M3 are chosen to resemble NSE for p = 
3.2 • 10 7 g cirr 3 , p = 5.0 • 10 7 g cnr 3 , andp = 1.0 • 10 8 g cm 4 . 
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Aq/10" erg g 4 
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Table 2. Parameters of the ashes for p = 3.2 • 10 7 g cm 



adding the difference of binding energies to the internal energy 
of the fuel (50% C, 50% O) and enforcing pressure equilibrium 
between the ashes and the enviroment. The latter condition is a 
consequence of the subsonic burning velocity of the deflagra- 
tion front. It guarantees the stationarity of the obstacles. 

Because the EOS in FLASH does not allow the use of inter- 
nal energy and pressure as state variables, we had to convert the 
conditions above into initial values for temperature and density. 
These are summarized in Tables|2] -|4] 



3.2. One-dimensional simulations 

The initial and boundary conditions for our analyzed system 
are the same as described in Sect. [2] Aditionally, we add an ash 
obstacle with width w at the position x = 512 cm (see Fig. [3}. 
We examined the behavior of the detonation front in this setup 
for widths w = 4 . . . 256 cm, compositions N,S, and M for the 
ash obstacle, and densities p = 3.2 • 10 7 ... 1.0 • 10 8 g cm" 3 . 

As measurement of the detonation strength, we analyzed 
the peak value of the pressure associated with the shock. After 
the phase of ignition with high overpressure, the peak value 
of the pressure converges to a stable state, as can be seen in 
Fig. [4] A stable shock front develops ahead of the reaction 
zone. The peak pressure declines steeply, as soon as the shock 
front hits the obstacle at the position x = 512 cm. For fuel 
densities of p = 3.2 ■ 10 7 gcm 4 , p = 5.0 ■ 10 7 gcm 4 , and 
p = 1.0 ■ 10 8 g cm -3 we found a reduction of peak pressure by 
1.5 ■ 10 25 erg cm 4 , 2.5 • 10 25 erg cm 4 , and 5.5 ■ 10 25 erg cm 4 . 

The composition seems to have virtually no influence on 
the outcome. This is surprising because, as one can see in Fig. 
|5] the ash is not a inert medium. Owing to strong local varia- 
tions in pressure and temperature, its composition changes sig- 
nificantly during the passage of the shock front. Nevertheless, 
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Fig. 4. Peak value of the pressure associated with the shock 
front over position for p = 5.0 • 10 7 g crrT 3 . The grey shaded 
areas illustrate the width of the obstacle from w = 4 - 256 cm. 
The legend shows the scaled width w s of the obstacles. 



independent of the initial composition, the final composition is 
nearly identical for all types of ash. 

Therefore ashes of type S and M generate nearly identical 
results in the analysis of the peak pressure, but only the ash of 
type N leads to a marginally greater decrease of the peak pres- 
sure. However, for obstacles with the smallest examined width 
w = 4 cm, the decline in the pressure maximum is approxi- 
mately 30% smaller than in the case of a width w = 16 cm or 
higher. But this difference only seems to have consequences in 
the case of the lowest analysed density of p = 3.2 • 10 7 g crrT 3 , 
where the detonation seems to reignite behind obstacles with 
widths w = 4 cm and w = 8 cm. This is shown in Fig.|6| where 
one can see that, immediately behind the obstacle, the pressure 
maximum increases to approximately 1.0 ■ 10 25 erg cm 4 above 
the value ahead of the obstacle. Afterward, the peak pressure 
quickly converges back to the characteristic pressure of a stable 
detonation. In all the other analyzed cases, the peak pressure 
also rises abruptly behind the obstacle, but a stable detonation 
front is not generated. In fact, the value of the peak pressure 
steadily decreases, which means that the shock front continu- 
ously loses strength. The shock front behind the obstacle is too 
weak to ignite a new self-sustaining detonation. 

3.3. Two-dimensional simulations 

For the 2D simulations, we set up a system with the same 
length, l x = 1024cm, as for the ID simulations. We set the 
width of the system to l y - 128 cm and initalized the system 
with 8 top-level blocks in the x-direction. Because each top- 
level block consists of 8 x 8 cells (or grid points), this ensures 
that the resolution in the x- and (/-directions were equal. We 
used the same boundary conditions in the x-direction as for 
the ID setup and set periodic boundary conditions in the y- 
direction. As in the ID simulation, we created an ignition area 
between x — cm and x — 25.6 cm that extended across the en- 
tire indirection. The inital conditions for the ignition area and 
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(a) p = 3,2 • 10 7 g cm 3 , ash type N 




(b) p = 5,0 • 10 7 g cm" 3 , ash type S 
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(c) p = 1,0 • 10 s g cm" 3 , ash type M2 

Fig. 5. Temporal evolution of the mass fractions of the six most 
important isotopes within the obstacle for an eviromental den- 
sity ofp = 5.0T0 7 g cm 4 . The legend shows the mass fractions 
of the specific isotopes at time t = 1,0 • 10~ 6 s. 



the enviroment (composition, temperature, and velocity of the 
fluid) were the same as in the ID simulation. 

As obstacles, we defined two semicircular areas of ash with 
radius r, centered at the positions (x, y) = (512,0) cm and 
(x, y) = (512, 128) cm (equivalent to a full circle of ash owing 
to the periodic boundary conditions, Fig. [7). 

We followed the evolution of the system with p = 5.0 • 
10 7 g crrT 3 and an ash obstacle of type M2 for different radii 
r. Limited computational resources restricted these simulations 
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front over position for p = 3.2 • 10 7 g cm -3 . The grey-shaded 
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The legend shows the scaled width w s of the obstacles. 
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Fig. 8. Carbon mass fraction using logarithmic scaling. The dis- 
tance between the two semicircles is 16 cm. 
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Fig. 7. Initial conditions for the 2D simulations with p — 5.0 • 
10 7 g cm -3 and a cirular obstacle of ash with radius r = 64 cm. 

to an effective resolution of 0.4 cm. The results of the 2D sim- 
ulations, therefore, have to be considered as preliminary. 

Figure [8] shows the development of the carbon mass frac- 
tion for a detonation front hitting a circular obstacle of ash with 
radius r — 56 cm, i.e., the distance a between the two semicir- 



cles of ash is 16 cm. One can see that the detonation squeezes 
between the two semicircles and re-expands after passing the 
narrowest point of the system. However, a "shadow" of un- 
burned material remains behind the former obstacle of ash. 
Figure [8] also shows a mushroom- like structure produced by 
the R ichtmyer-Meshkov instability ( Richtm ved 1 960t [Meshkovl 
1970) occuring when a shock front passes the contact discon- 
tinuity between two fluids of different densities. Nethertheless, 
the detonation front is not stopped by an obstacle of these pro- 
portions. 

A different result was found for the distance a < 12 cm, 
as can be seen in Fig. |9] Although the detonation front prop- 
agates through the small funnel between the two semicircles 
of ash, it becomes unstable and dies right behind the obsta- 
cle. Apparently, the funnel of unburned material between the 
two obstacles is too small to provide enough energy for a self- 
sustaining detonation front. 

4. Conclusions 

Our results have immediate consequences for models of ther- 
monuclear detonations in large-scale SN Ia simulations. In ad- 
dition to an accurate reproduction of the detonation speed, 
which is usually achieved by simulating an unresolved det- 
onation that is directly coupled to the nuclear reactions via 
the coarse-grained temperature, we must enforce the new re- 
quirement that the front stops immediately when it encoun- 
ters burned material, since the typical size of ash regions that 
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it may cross can never be resolved in practice. Whether or 
not unresolved, temperature-coupled detonations (frequently 
used in multi-dimensional explosion models on large scales) 
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and ignition. If the spatial resolution is chosen too coarsely in 
Eulerian methods, numerical diffusion of heat from the ash pre- 
conditions the fuel and helps to re-ignite the detonation (ID 
Lagrangian codes are not subject to this restriction). We have 
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a numerical artifact under closer scrutiny. 
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inhomogeneities over a wide range of length scales, whereas 
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ison of 3D simulations of pure deflagration and delayed det- 
onation models for SNe la will reveal the ability of delayed 
detonations to enhance the chemical stratification of the ejecta 
and to reduce the amount of unburnt material at low velocities. 
This investigation is currently in progress. 



Colella, P. & Glaz, H. M. 1985, Journal of Computational Physics, 59, 
264 

Colella, P. & Woodward, P. R. 1984, Journal of Computational 

Physics, 54, 174 
Fryxell, B., Olson, K„ Ricker, P., et al. 2000, ApJS, 131, 273 
Gamezo, V. N, Khokhlov, A. M., & Oran, E. S. 2005, Astrophys. J., 

623, 337 

Golombek, I. & Niemeyer, J. C. 2005, A&A, 438, 61 1 

Khokhlov, A. M., Oran, E. S., & Wheeler, J. C. 1997, ApJ, 478, 678 

Kozma, C. et al. 2005 

Leonard, D. C, Li, W., Filippenko, A. V., Foley, R. J., & Chornock, 
R. 2005 

Livne, E. 2004, in Cosmic explosions in three dimensions, 110 
Mazzali, P. A. et al. 2005, Astrophys. J., 623, L37 
Meshkov, E. 1970, NASA Tech. Trans., F13, 074 
Niemeyer, J. C. 1999, ApJ, 523, L57 
Niemeyer, J. C. & Woosley, S. E. 1997, ApJ, 475, 740 
Richtmyer, R. 1960), Commun. Pure Appl. Math., 13, 297 
Riess, A. G. et al. 2005 

Ropke, F. K., Gieseler, M., Reinecke, M., Travaglio, C, & Hillebrandt, 
W. 2005 

Sharpe, G. J. 1999, MNRAS, 310, 1039 
Timmes, F. X. 1999, ApJS, 124, 241 

Timmes, F. X., Hoffman, R. D., & Woosley, S. E. 2000a, ApJS, 129, 
377 

Timmes, F. X. & Swesty, F. D. 2000, ApJS, 126, 501 
Timmes, F. X., Zingale, M., Olson, K., et al. 2000b, ApJ, 543, 938 
Zingale, M., Woosley, S. E., Rendleman, C. A., Day, M. S., & Bell, 
J. B. 2005 



Acknowledgements. We would like to thank Marcus Bruggen and 
Pawel Ciecielag for their help with the FLASH code, and Christian 
Klingenberg, Irina Golombek, Ewald Miiller, Jan Pfannes, and 



